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ABSTRACT 

We explore the evolution of the emissions by accelerated electrons in shocked shells driven by jets 
in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, we 
evaluated the broadband emission spectra by properly taking into account adiabatic and radiative 
cooling effects on the electron distribution. The synchrotron radiation and inverse Compton (IC) 
scattering of various photons that are mainly produced in the accretion disc and dusty torus are 
considered as radiation processes. We show that the resultant radiation is dominated by the IC 
emission for compact sources (< lOkpc), whereas the synchrotron radiation is more important for 
larger sources. We also compare the shell emissions with those expected from the lobe under the 
assumption that a fractions of the energy deposited in the shell and lobe carried by the non-thermal 
electrons are e e ~ 0.01 and e e! i bc ~ 1, respectively. Then, we find that the shell emissions are 
brighter than the lobe ones at infra-red and optical bands when the source size is > lOkpc, and the 
IC emissions from the shell at > 10 GeV can be observed with the absence of contamination from 
the lobe irrespective of the source size. In particular, it is predicted that, for most powerful nearby 
sources (L- 3 ~ 10 47 ergs s _1 ), ~ TeV gamma-rays produced via the IC emissions can be detected by 
the modern Cherenkov telescopes such as MAGIC, HESS and VERITAS. 

Subject headings: particle acceleration — radiation mechanisms: non-thermal — galaxies: active — 
galaxies: jets — 



1. INTRODUCTION 

It is well established that radio-loud active galactic nu- 
clei (AGNs) are acco mpanied by relativistic jets (e.g., 
iBegelman et al.lll984 for review). These jets dissipate 
their kinetic energy via interactions with surrounding in- 
terstellar medium (ISM) or intracluster medium (ICM), 
and inflate a bubble composed of decelerated jet matter, 
which is often referred to as cocoon. Initially, the cocoon 
is highly overpressured a gainst the ambient ISM/ICM 
(jBegelman fc Cio"m1ll989| ) and a strong shock is driven 
into the ambient matter. Then a thin shell is formed 
around the cocoon by the compressed ambient medium. 
The thin shell structure persists until the cocoon pres- 
sure decr eases and the pressure equilibrium is eventually 
achieved ([Reynolds et alJfeOOll ). 

While large number of radio observations identified co- 
coons with the extended radio lobe, no clear evidence 
of radio emissions is found for the shocked shells (e.g., 
iCarilli et al.l I1988T) . Due to the lack of detections, in 
the previous studies on the extragalactic radio sources, 
it is usually assumed that non-thermal emissions are 
dominated by or o riginated only in the cocoon (e.g., 
IStawarz et al.l 120081) . However, since strong shocks are 
driven into tenuous ambient gas with a high Mach num- 
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ber, the shocked shells are expected to offer site of par- 
ticle acceleration as in the shocks of supernova remnants 
(SNRs) and theref ore give ri se non-thermal emissions 
(|Fuiita et all 120071 : iBerezhkol I2008f ). Hence, although 
the observations at radio seems to be unsuccessful, non- 
thermal emissions from the shell may be accessible at 
higher frequencies. In fact, while not detected in radio, 
recent deep X-ray observation have reported the presence 
of non-thermal emissions from the shell associated with 
the radio galaxy Centaurus A ([Croston et al. 200j|). 

Theoretically, emissions by the accelerated particles re- 
siding in shocked shells have been studied bv lFuiita et al.1 
(2007). However, they paid attention only to the the ex- 
tended sources of ~ 100 kpc and the inverse Compton 
(IC) scattering of external photons was not included in 
the radiative processes. The cooling effects on the energy 
distribution of non-thermal electrons were not consid- 
ered, either. Motivated by these backgrounds, we explore 
in this paper the temporal evolution of the non-thermal 
emissions by the accelerated electrons in the shocked 
shells, properly taking into account the Comptonization 
of photons of various origins as well as the cooling effects 
on the electron distribution. Focusing on the powerful 
sources which host luminous quasar in its core, we show 
that the shell can produce prominent emissions ranging 
from radio up to ~ 10 TeV gamma-ray and discuss the 
possibility for the detection. 

The paper is organized as follows. In Sj2j we introduce 
the dynamical model, which describes the evolution of 
the shell expansion, and explain how the energy distri- 
bution of electrons residing in the shell and the spectra of 
the radiations they produce are evaluated based on the 
dynamical model. The obtained results are presented in 
S3 In 21 we compare the emissions from the shell with 
those from the lobe and discuss its detectability. We 
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close the paper with the summary in <j5j 

2. MODEL FOR EMISSIONS FROM THE SHELL 
2.1. Dynamics 

In this section, we give the model to approximately de- 
scribe the dynamics of the expanding cocoon and shell 
and to provide a basis, on which the energy distribu- 
tion of accelerated electrons and the radiation spectra 
are estimated. The schematic picture of the model is 
illustrated in Fig. [TJ For simplicity we neglect the elon- 
gation of cocoon and shell in the jet direction and assume 
that they are spherical. We also assume that the shell 
width, SR, is thin compared with the size of the whole 
system, R, and ignore the difference between the radii of 
the bow shock and the contact surface. This assumption 
is valid as long as the expansion velocity ha s a high Mach 
number (see e.g.. lOstriker fc McKeel Il988l ) . We further 
assume that the kinetic power of jet, Lj, is constant in 
time on the time scale of relevance in this study. 

Under the above assumptions, the dynamics of the ex- 
panding cocoon and shell can be appr oximately describe d 
by the model of stellar wind bubbles (Castor c t al.lll975l ). 
The basic equations are the followings: ([T]) the equation 
for the momentum of the swept-up matter in the shell 
and (0) the equation for the energy in the cocoon. They 
are expressed, respectively, as 

j t (M s (t)R(t)) = 4irR(t) 2 P c (t), (1) 



d_ fP c (t)V c (t) 
dt 



7c 



- 1 



Pc(t) 



dVcjt) 
dt 



(2) 



where R = dR/dt, V c = 4irR 3 /3, and P c are the expan- 
sion velocity of the shell, the volume and pressure of the 
cocoon, respectively, and % is the specific heat ratio for 
the plasma inside the cocoon. The swept-up mass in the 
shell is defined as M s = J fl Anr 2 p a {r)dr with the mass 
density of the ambient medium, p a . 

In this study, p a is assumed to be given by p a (r) — 
Po{r/ro)~ a with ro and po being the reference position 
and the mass density there, respectively. Then Eqs. (TTJ 
and ([2]) can be solved analytically and the solution is 
given as 



R(t) = C R r 



a/ (a— 5) 



l/(5-a) 



t 3/(5-a) j (3) 



where Cr is a function of a and 7 C given as 



(3-a)(5-q) 3 (7 c ~l) 
47r{2a 2 + (1 - 187 c )a + 63f G - 28} 



l/(5-o) 



Employing the non-relativistic Ran kine-Hugoniot con- 
dition in the strong shock limit ([Landau fc Lifshit3 
1959), which is justified by the assumption that the shell 
expands with a high Mach number, we obtain from Eq. 
(|3j) the density, p s , and pressure, P Sl in the shell, which 
are assumed to be uniform, as follows: 
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Fig. 1. — The schematic picture of the model employed in this 
study to approximate the expansion of the cocoon and shell pro- 
duced by the jet pushing into the ambient medium from AGN. 



where 7 a is the specific heat ratio for the ambient 
medium. By equating M s with the mass in the shell 
given by V s p s , where V s = 4irR 2 SR is the volume of the 
shell, the shell width is obtained as SR = (% — + 
1)(3 — a)]~ 1 R . It is found that the ratio of SR to R 
does not depend on time and the thin-shell approxima- 
tion is reasonably good for the typical values, 7 a = 5/3 
and < a < 2. 

The total internal energy stored in the shell, E s = 
PsVs/ila — 1), is one of the most important quantities, 
since it determines the energy budget for the radiation. 
From Eqs. © and E s can be expressed as 



E s - /£jt, 



(5) 



where / is a fraction of the total energy released by the 
jet (Ljt), which is converted to the internal energy of the 
shell, and is given by 



/ = 



18( 7c -l)(5-a) 



{% + l) 2 [2a 2 + (1 - 18 7c )a + 63 7c - 28] 



(6) 



which turns out to be time-independent. For typical 
numbers, j c = 4/3, % — 5/3, < a < 2, / depends 
on a only weakly and / ~ 0.1. 

2.2. Energy Distribution of Electrons 

Assuming that a fraction, e e , of the energy deposited 
in the shell, E s , goes to the non-thermal electrons, we 
evaluate their energy distribution. We solve the kinetic 
equation governing the temporal evolution of the energy 
distribution of electrons, N{^ e , t), which is given as 



dN{ lei t) 
dt 



a 



( 7e )iV(7e,*)]+Q(7e 



(7) 



where j e , 7 CO oi(7e) = —dj e /dt, and Q(7 e ) are the Lorentz 
factor, the cooling rate via adiabatic expansions and ra- 
diative losses, and the injection rate of non-thermal elec- 
trons, respectively. The latter two, that is, the injection 
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rate Q{^ e ) and the cooling rate 7 CO oi, which will be de- 
scribed in detail below, are evaluated based on the dy- 
namical model described in the previous section. 

Following the theory of the diffusive shock accelera- 
tion (DSA) (IBelllll978l : [Dmr^[T98llBlandford fc Eichlerl 
fl987h . we assume that the non-thermal electrons are in- 
jected into the post-shock region with a power-law energy 
distribution of the form 



Q(7e) = Kj e p for 7 min < 7 e < 7,, 



(8) 



where 7 m i n and 7 max correspond to the minimum and 
maximum Lorentz factors, respectively. In this study 
we set 7mi n = 1 and the power-law index, p, is fixed 
to 2, an appropriate value for the linear diffusive shock 
acceleration in the strong shock limit. The maximum 
Lorentz factor is obtained by equating the the cooling 
rate, 7 CO oi: to the acceleration rate given by 



7a 



3eBR 2 
20£m e c 3 ' 



(9) 



where B is the magnetic field strength in the shell. The 
so-called "gyro-factor" , £, can be identified with the ra- 
tio of the energy in ordered magnetic fields to that in 
turbulent ones (£ = 1 is the Bohm limit). The nor- 
malisation factor, K, in the injection rate is determined 
from the assumption that a fraction, e e , of the shock- 
dissipated energy is carried by the non-thermal electrons: 
X^"(7e ~ l)m e c 2 Q(-/ e )d-y e = e e dE s /dt = fe e L y A 
rough estimation of the normalization factor is obtained 
as K ~ 0.1e e Lj/[m e c 2 ln(7 max )], where we used / ~ 0.1 
( §2.1|) . It is noted that since the factor K is proportion- 
ate to e e and Lj , the resultant luminosity of non-thermal 
emissions also scales in the same manner with these quan- 
tities. 

In the cooling rate, 7 CO ol> both radiative and adiabatic 
losses are taken into account. As for the former, the 
synchrotron radiation and the IC scattering off photons 
of various origins. The cooling rate for the adiabatic 
dynamical expansion is expressed for electrons with a 
Lorentz factor 7 e as 



7ad 



1 R 

3R 



7e 



(10) 



The rates of the radiative coolings via the synchrotron 
radiation and IC emissions are given, respectively, as 



7syn 



A(J T U B 2 

3m„c 



and 



7IC = -5 7e-FKN(7e), 

6m e c 



(11) 



(12) 



where ax, c and m e are the cross section of Thomson 
scattering, the speed of light and the electron mass, re- 
spectively. The energy densities of magnetic fields and 
photons in the shell are denoted by Ub = B 2 /8tt and 
Uph, respectively. The function i 7 KN(7e) encodes both 
the distributions of seed photons and the Klein-Nishina 
(KN) cross section for the Compton scattering and re- 
duces to unity in the Thompson limit. Note, however, 
we do not employ this limit but calculate F^ile) from 



the differential cross-section for IC scattering given in 
IBlumenthal fc Gouldl <fl970h . 

The typical magnetic field strengths in el- 
liptical galaxies and cl usters of galax ies are 
a few mG (e.g.. IMoss fc Shukurovl 119961: 



Vikhhnin. Markevitch fc Murray! 120011: iClarke et alJ 



200H ICarilli fc Tavlorl 12002b ISchekochihin et all I200HK 
Assuming that as it passes through a shock wave, the 
magnetic field is adiabatically compressed by a factor of 
1 ~ 4 depending on the obliqueness of the shock to the 
magnetic field lines, we choose B — 10/iG as a fiducial 
value for the magnetic field strength in the shell when 
evaluating the acceleration rate, 7 aC cci, and synchrotron 
cooling rate, 7 syn . The corresponding energy density of 
magnetic field is given by 



U B « 4.0 x lQ- r2 B 2 _ 5 erg cm" 



(13) 



where B_ 5 = B/10fiG. 

In evaluating the cooling rate for IC scattering, 71c, 
we take into account various seed photons of re l evance 
in this context. In their paper, iStawarz et all J2008) 
explored high energy emissions by relativistic electrons 
in the radio lobes through the IC scatterings of various 
photons. Among them are UV emissions from the accre- 
tion disc, IR emissions from the dusty torus, stellar emis- 
sions in NIR from the host galaxy and synchrotron emis- 
sions from the radio lobe. In addition to these emissions, 
we also consider CMB as seed photo ns in the present 
study. Following IStawarz et al.l (|2008l ). we assume that 
the photons from the disc, torus, host galaxy, and CMB 
are monochromatic and have the following single fre- 
quencies: v\jy = 2.4 x 10 15 Hz, v m = 1.0 x 10 13 Hz, 
fNiR = 1.0 x 10 14 Hz, and v CMB = 1.6 x 10 11 Hz. The 
photons from the radio lobe are assumed to have a con- 
tinuous spectrum given by L„ \ or)e cx i/~ a75 . 

From the luminosity, L\jy, of the UV emissions by the 
accretion disc, the energy density of these photons in the 
shell is given approximately as 



0uv = 



L 



uv 



4irR 2 c 
3.0 x 10" 



^uv,46#i 2 erg cm" 



(14) 



where Luv,46 = ^uv/10 46 erg s \ an d R\ = i?/lkpc. 
As is assumed in IStawarz et al.1 (|2008[) , we take L uv = 
10 46 ergs s _1 for sources with jet kinetic power of Lj > 
10 45 ergss"" 1 , and Luv = 10 45 ergs s _1 for Lj < 
10 45 ergs s" 1 as fiducial values. The adopted values are 
appropriate for sources hosting a luminous quasar. Ac- 
cording to the broadband observa tions of quasars (e.g., 
lElvis et alJH99l Uiang et al.|[2006l) , the IR and UV emis- 
sions are approximately comparable although there are 
some variations in the relative strengths from source to 
source. We hence assume that the luminosity of the IR 
emissions from the torus, £ir, is equal to that of the UV 
emissions from the disc, L\jy, and the energy density of 
IR photons in the former case is estimated as 



U m 



L 



IB 



3.0 x 10 _9 L 



Hi 46-Ri 2 



erg cm 3 , (15) 



4itR 2 c 

where Lir,46 = -^ir/10 46 erg s~ 4 . In evaluating the en- 
ergy density of photons from the host galaxy, we assume 
that most of the photons are emitted by stars in the the 
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core region with the radius of ~ 1 kpc (|de Ruiter et al.l 
2005). Considering only the sources with R > kpc, we 
can estimate the energy density of the optical photons 
from the luminosity, Lnir, as 



t r NIR = 



AttR 2 c 

3.0 x 10~ 10 L NIFM5 i?r 2 erg cm' 



(16) 



where Lnir,45 = £Nm/10 45 ergs s 1 . The energy density 
of CMB photons is given by 



Ucmb « 4.2 x 10 13 erg cm 



(17) 



The redshift is ignored for simplicity. Finally, the energy 
density of photons emitted from the lobe is obtained by 
assuming that a fraction rj of the jet power is radiated as 
radio emissions from the lobe (i.e., J Ly^^dv = rjLj). 
Assuming rj — 10 -2 as a fiducial case, it is given as 

Ulobc ~ 4nR2c 

«3.0 x 10- 12 T]_ 2 L i5 R^ 2 erg cm' 3 , (18) 

where r/_ 2 = i]/10~ 2 and L 45 = Lj/10 45 ergs s _1 . 

The cooling rate for the IC scattering is obtained from 
Eq. ([T2j) by plugging in the energy densities U p h ob- 
tained above and calculating i*kN(7e) f° r the- assumed 
photon distributions. Each contribution from the disc, 
torus, host galaxy, CMB, and lobe is denoted as 71c, uv, 
Tic.ir, 7ic,nir, 7ic,cmb, and 7ic,lobe- The total cooling 
rate 7 CO oi is the sum of the rates for the adiabatic, syn- 
chrotron and IC losses. (7 coo i = 7ad + 7syn + 7ic,uv + 
7ic,nir + 7ic,cmb + 7iC,lobe)- Note that the synchrotron 
self-Compton is ignored, since its effect is negligible in 
any case considered in the present study. 

Now that the injection rate, Q("f e ), and the cooling 
rate, 7 CO oi, have been evaluated, the energy distribution 
of non-thermal electrons, N{^ e , t), is obtained by putting 
these quantities in Eq. ([7]). Although Q(j e ) is time- 
dependent through 7 max and so is 7 CO oi because of the 
power-dependence of the expansion rate (see Eq. d3J)), we 
ignore these variations over the dynamical time scale ~ t. 
Then, employing the instantaneous values of 7 CO oi(7e) 
and Q(7e) evaluated at each time and fixin g them, we can 
solve Eq. Q by the L aplace transforms (jMelrosd 119801 : 
iManolakou et al.l 12007ft and the solution can be written 



JV(7e,*) 



1 



7u 



dl' e Q(l' e ), (19) 



7cool(7e) J ls 

where the upper limit of the integral is given by 



7u 



7*(7e 
7max 



fo 



r 7n 



< 



7< 



< 



for 7b r < 7 e < 7 1J 



7br, 



Here 7* and 7t> r are determined by the relation 
r (7*,7e) = T(7 max ,7b r ) = t, where the function 

r (7'i7) = /~ ^7" licoo\{l") gives the time it takes an 
electron to cool from 7' to 7 in the Lorentz factor. Note 
that 7br corresponds to the Lorentz factor, where the en- 
ergy distribution of non-thermal electrons shows a break 
owing to the radiative coolings (see ^3.1l for detail). 



2.3. Radiation spectra 

From the energy distribution of non-thermal electrons 
just obtained, the spectrum of the synchrotron radiation 
is calculated as 



-P S yn(f,7e)Af(7 e )<i7e 



(20) 



Here, P syn (^, 7 e ) is the pitch-angle-averaged power spec- 
trum for a single electron given by 



-P S yn(^,7e) 



V3e 3 B 
2m e c 2 



sinW I — ) d0, 



where is the pitch angle and v c = 3ei?7 2 sin0/47rm e c 
is the characteristic frequency of the emitted photons. 
The function F(x) is defined by F(x) = x J^ K 5 / 3 (y)dy, 
where K 5 / 3 (y) is the modifie d Bessel function of the 5/3 
order (Rybick i fc Lightmanl 119791) . Note that the syn- 
chrotron self-absorption is ignored, since it is important 
only at low frequencies below 10 7 Hz in the current con- 
text. 

In the same way, the spectrum of the IC scattering is 
obtained as 



L 



u,IC 



Plc(*7 7e)iV(7e)(27e 



(21) 



where 



Pic(v,7e) = hvc I n ph (u s )aic(iy,iy a ,^ e )dv s . 
Jo 

is the power spectrum for a single electron in isotropic 
photon fields. Here, n p h(i / s) and 010(27 v B , 7e) are the 
number density of seed photons per unit frequency and 
the differential cross s e ction for the IC scattering given by 
iBlumenthal fc Gould! (|1970| ). which is valid both in the 
Thompson and KN regimes. The origins of the seed pho- 
tons have been discussed in §2.21 The number densities 
per unit frequency of these photons are then given as fol- 
lows: n p h(^ s ) = (Uph/hVph)S(is s — Vph) for the UV emis- 
sions from the accretion disc (J7 p h = UuVyVph = v\jv), 
the IR emissions from the torus (C/ p h = C^m, v p \i = 
vm), the NIR emissions from the host galaxy (U p h — 
Uv,v p h = Vv), and CMB (U ph = U C MB,v P h = ^cmb); 
^ph(i / s) = Lv^obe/liVsAir R 2 c for the lobe emissions. 

3. BROADBAND SPECTRUM OF THE SHELL 

In the following we focus on powerful radio sources 
(Lj ~ 10 45 — 10 47 ergs s -1 ) hosting luminous quasars and 
present the temporal evolutions of the energy distribu- 
tion of non-thermal electrons and the resultant radiation 
spectra. 

As a fiducial case, we set the parameters for the am- 
bient matter as po = 0.1m p cm~ 3 ,ro = 1 kpc and 
a = 1.5, i.e., p a {r) — 0.1m p po.i( r /lkpc) _1 ' 5 cm -3 , where 
m p is the proton mass and po.i = Po/0.lm p . These 
are indeed typical values for elliptical galaxies (e.g. , 
Mulchaev fc Zabmdofflll998tlMathews & Brighentill2003l: 
Fukazawa. Makishima. fc Ohashil 120041 : iFukazawa et al.1 
2006J). We adopt % = 4/3 and % = 5/3. Then, R is 
given from Eq. © as 



R(t) 



„„ -2/7 r 2/7,6/7 , 
22 Po.l L 45 W k P c > 



(22) 
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where tj = t/10 7 yr. As mentioned in $2.21 the power- 
law index, p, for the injected non-thermal electrons (see 
Eq. (JSJ) is fixed to 2. Although the parameters £ and 
e e which characterize the electron acceleration efficien- 
cies are highly uncertain, it is natural to expect that 
the range of these values are similar to those observed 
in the SNRs, since the nature of shocks considered here 
resembles those in the SNRs in that they are strong non- 
relativistic collision-less shocks driven into ISM. Here 
we take £ = 1, which corresponds to the Bohm diffu- 
sion limit, and e e = 0.0 1 as illustrat ive values, based 
on the the observations of SNRs (e.g., iDver et al.ll2001t 
Ellison et all 120011: iBamba et al.l 120031 lYamazaki et al.l 
20041: iStage et al.l 120061: iTanaka et al.l 12008ft . The lumi- 
nosities of the emissions from accretion disc and dust 
torus are taken as Lyv = -^ir = 10 46 ergs s _1 for 
Lj > 10 45 ergs s _1 , and i uv = Im = 10 45 ergs s _1 
for Lj < 10 45 ergs s _1 , while fixed value of Lnir = 
10 45 ergs s _1 is adopted for the emissions from the host 
galaxy. The parameters for the magnetic field and lobe 
emissions are chosen as B = 10/iG and r\ = 0.01. The 
employed values of the parameters are summarized in 
Table [TJ 

3.1. Evolution of Non-thermal Electrons 

In Figs. [2] and El we display the cooling and ac- 
celeration time scales (left panels) and the energy dis- 
tribution of non-thermal electrons (right panels). The 
jet powers are chosen to be Lj = 10 45 ergs s _1 and 
Lj = 10 47 ergs s _1 in Figs. [2]and[3j respectively. The 
top, middle and bottom panels in the figures are given for 
the source sizes of R = 1, 10, and 100 kpc, respectively, 
which in turn corresponds to the different source ages. 
In addition to the total cooling time scale, t coo \, we show 
the contributions from the adiabatic loss, t a( j = 7e/7ad, 
synchrotron radiation, t syn — r y e /isyni an d IC scatter- 
ings of the UV disc photons, iic,uv = 7e/7ic,uv, IR 
torus photons, iic,iR = 7e/7ic,iR, NIR host galaxy pho- 
tons, iic,NiR = 7e/7ic,NiR, CMB photons, £i C ,cmb = 
7e/7ic,CMB, and lobe photons, iicjobc = 7e/71C,lobe, to- 
gether with the source age, t. 

As mentioned in W2.2\ the Lorentz factor at the spec- 
tral break, 7b r , corresponds to the electron energy, above 
which cooling effects become important. It is hence de- 
termined roughly by the condition t ~ t coo \ as shown 
in the left panels of Figs. [2] and [3] Since the adiabatic 
loss is dominant over the radiative coolings only below 
7br, the latter is always more important when the cool- 
ings affect the electron distribution. The relative im- 
portance of the synchrotron emissions and IC emissions 
depends on the electron energy as well as the energy den- 
sities of magnetic fields and photons. When the source is 
young and hence small, the energy loss is dominated by 
the IC emissions, since the energy density of photons is 
larger than that of magnetic fields (referred to as the IC- 
dominated stage). As the the source becomes larger, on 
the other hand, the energy density of photons decreases 
(C/ph oc R~ 2 ) and the synchrotron loss becomes more im- 
portant (the synchrotron-dominated stage). In the initial 
IC-dominated stage, the break Lorentz factor increases 
with the source size as 7b r <x i? 5 / 6 , where we used the 
relations t oc i? 7 / 6 and t coo \ oc 7 c T 1 ?/ ph 1 oc 7 ( T 1 i? 2 . On 
the other hand, 7b r decreases with the source size as 



7br oc R 7 / 6 in the synchrotron-dominated stage, since 
the relation t coo i oc 7 c T 1 f/B 1 oc 7,T 1 i? holds. Regard- 
ing the dependence on the jet power, for a given source 

1/3 

size R, 7br increases with Lj roughly as 7b r oc Lj be- 
cause the dependencies of the source age and the radia- 
tive cooling time scale on the jet power can be expressed 

as t oc *' 3 and 7b r oc 7 ( T 1 i!\ respectively. 

Among the contributions to the IC losses, the scat- 
tering of the IR torus photons is the largest at least in 
the IC-dominated stage thanks to the high energy den- 
sity of the IR photons. Although the UV disc photons 
are assumed to have the same energy density as the IR 
photons (Utjv = the cooling by the former is sup- 

pressed for 7 e > m e c 2 / Ahv\j\j ~ 1.3 x 10 4 by the KN 
effect as seen in Figs. [2] and [3J As a result, unless Luv 
exceeds Lir significantly, most electrons with 7 e > 10 4 
cool predominantly through the IC scattering of IR pho- 
tons. Then the transition from the IC-dominated stage 
to the synchrotron-dominated stage occurs roughly at 

i?ic/ sy n ~ 27L I 1 / 2 46 BZ 5 2 kpc, (23) 

which corresponds to the condition J7ir ~ Ub- Since the 
IC scattering of IR torus photons is also suppressed by 
the KN effect above the Lorentz factor given by 

7kn = m e c 2 /Ahv m - 3 x 10 6 , 

the coolings by the synchrotron emissions and/or the IC 
scattering of lobe photons are important at the high en- 
ergy end of non-thermal electrons (7kn < 7e < ) for 
the source with R < -Ric/syn (see, e.g., the middle left 
panel of Fig. [2] and the top and middle left panels of Fig. 
[3]). For even larger sources with R > -Ric/sym the syn- 
chrotron emissions are the dominant cooling mechanism 
at all energies. The contributions from the IC scatterings 
of CMB and host galaxy photons, on the other hand, are 
modest at most during the whole evolution. 

The maximum Lorentz factor of the injected electrons, 
which is determined by the condition t coo i — t accc \ (or 
equivalently 7 coo i = 7accd), is 7max ~ 10 7 - 10 9 for 
the source sizes R ~ 1 — 100 kpc and the jet power 
Lj ~ 10 45 — 10 47 ergs s _1 . Regarding the dependence 
on the source size, 7 max increases with R for compact 
sources in which the cooling at 7 e ~ 7max is dominated by 
the IC emissions. This is because the cooling time scale 
increases rapidly with size (i coo l oc 7 f r 1 LkN(7e) _1 C^ p i 1 1 oc 
7~ 1 i 7 KN(7e)~ 1 -R 2 ), in contrast with the slow increase in 
acceleration time scale (i a ccei oc "/ e R~ 2 oc 7 e i? 1 / 3 ). On 
the other hand, for larger sources in which the cool- 
ing at 7 e ~ 7max is dominated by the synchrotron 
emissions, 7 max decreases slowly with the source size 
(7max oc i? -1 / 6 ), since the cooling time scale is indepen- 
dent of the size (t coo i oc 7 < T 1 C/g 1 oc 7 ( T 1 i?°). Note that 
this transition takes place earlier than the transition from 
the IC-dominated stage to the synchrotron-dominated 
stage mentioned above owing to the suppression of IC 
emissions in the high energy electrons (-f e > 7kn) by the 
KN effect. For a given source size, higher values of j max 
are obtained for sources with larger jet powers, since the 

acceleration time scale is shorter (£ a ccci oc 7 e Lj 2//3 ). 

The resultant energy distributions of non-thermal elec- 
trons, iV(7 e )j are understood as follows. They have a 
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TABLE 1 
List of the Model Parameters. 



Parameters Symbols Employed Values 



power-law index for density profile of ambient matter, p a (r) 


a 






1.5 


reference position in p a (r) 


ro 






1 kpc 


mass density, p a (ro), at r = ro 


PO 






0.1m p cm -3 


magnetic field strength 


13 






10 fJ,G 


luminosity of UV emissions from accretion disc 


iuV 


10 46 erf 


;s s~ 


1 (for Lj > 10 45 






10 45 erf 


;s s~ 


1 (for Lj < 10 45 


luminosity of IR, emissions from dust torus 




10 46 erg 


;s s~ 


1 (for Lj > 10 45 






10 45 erf 


;s s~ 


1 (for Lj < 10 45 


luminosity of NIR emissions from host galaxy 


^NIR 






10 45 ergs s _1 


ratio of lobe luminosity to jet power 


V 






0.01 


power-law index for energy distribution of injected electrons 


V 






2 


gyro-factor 


£ 






1 


energy fraction of non-thermal electrons 








0.01 



ergs s ) 
ergs s ) 



sharp cut-off at j e = 7 max whereas the original spec- 
trum N(j e ) oc 7~ 2 of the injected electrons is retained 
at low energies below 7b r - Since the cooling effects 
are negligible, the energy spectrum for j e < 7b r are 
roughly given by N(j e ) ~ Q(j e )t. The spectrum steep- 
ens for 7 e > 7b r owing to the radiative cooling and is 
roughly determined by the energy dependence of t coo \ as 
N(le) ~ Q(7e)^cooi(7e)- As mentioned above, for sources 
smaller than -Ric/sym the IC scattering of IR torus pho- 
tons determines the spectral shape up to ~ 7kn, where 
the KN effect kicks in and hardens the energy spectra for 
higher energies (see, e.g., the top and middle panels of 
Figs. [2] and El). This feature is absent for larger sources 
(R > i?ic/ S yn)j because the synchrotron loss (j S y n oc j 2 ) 
is dominant at all energies and the power is simply re- 
duced by 1 (N(j e ) oc 7jT 3 ) from the original spectrum 
of the injected electrons. For given source size, higher 
values of N(j e ) are obtained for sources with larger Lj, 
since the electron injection rate Q(j e ) is higher. 

3.2. Evolution of Radiation Spectra 

In Fig. U]we show the radiation spectra, vL v , which 
are obtained by inserting into Eqs. (l2Cfl) and (l2Tj) the 
energy distributions of non-thermal electrons in the pre- 
vious section. It is assumed that the jet powers are 
Lj = 10 45 ergs s _1 (left panels) and Lj = 10 47 ergs s _1 
(right panels). The top, middle and bottom panels of the 
figure correspond to the source sizes of R = 1 kpc, 10 kpc 
and 100 kpc, respectively. In addition to the total lumi- 
nosity (thick solid line), we show the contributions from 
the synchrotron emissions (thin solid line) and the IC 
scatterings of UV disc photons (long-short-dashed line), 
IR torus photons (dot- dashed line), NIR host galaxy pho- 
tons (dotted line), CMB photons (long-dashed line) and 
lobe photons (short-dashed line). 

The synchrotron emissions are the main low-frequency 
component, which extends from radio to X-ray ~ 
2(7max/10 8 ) 2 -B_ 5 keV. The IC emissions become re- 
markable at higher frequencies up to gamma-ray ~ 
50(7max/10 8 ) TeV. As mentioned in ^3.11 the dominant 
radiative process changes from the IC emissions to the 
synchrotron radiations as the source becomes larger. As 
a result, the former is more luminous than the latter 
when the source is young and compact (R < -Ric/syn) and 
vice versa. Owning to the hardening of the energy distri- 
butions of non-thermal electrons at -f e > 7kn (see Figs. 
[2] and [3]) for compact sources (R < L? IC / syn ), their syn- 
chrotron spectrum becomes also harder than for larger 



sources. This feature is not so remarkable for most of 
the IC components mainly because the emissions them- 
selves are suppressed by the KN effect. The IC scattering 
of IR dust-torus photons dominates over other IC com- 

1/2 

ponents up to the source size of R ~ 85L T p, 46 kpc. For 
larger sources, on the other hand, the IC scattering of 
CMB photons becomes more important, since CMB has 
the largest energy density of all photons considered in 
this paper. The contributions from the UV disc photons 
and NIR host galaxy photons are modest at best through 
the entire evolution. 

The energy injection rate above 7b r is roughly equal 
to the radiative output ( r y 2 m e c 2 Q("f e ) ~ vL v ) of non- 
thermal electrons in this energy regime because the cool- 
ing time scale is shorter than the dynamical time scale. 
Since the energy injection rate is independent of the elec- 
tron energy (^ 2 Q(^ e ) oc 7°), these non-thermal electrons 
produce a rather flat and broad spectrum (vL u oc v°) 
in the corresponding frequency range. From the relation 
7 2 m e c 2 Q(7 e ) = Km e c 2 ~ 0.1e e Lj/[ln( 7max )] (see gOJ), 
we can give a rough estimate to the peak luminosity as 

(vL v ) veak ~ 4.0 x 10 40 e_ 2 i45 ergs s" 1 , (24) 

where e_2 = e e /0.01. In the above equations, we ignored 
the dependence of luminosity on 7 max because the lu- 
minosity only scale logarithmically with its value, and 
employed a typical value 7 max ~ 10 8 . The feature is 
clearly seen in Fig. |4] Indeed, the spectra are flat with 
a peak luminosity given approximately by Eq. (|24j) . As 
mentioned in i)2.21 the emission luminosity scale approx- 
imately linearly with the acceleration efficiency e e and 
the jet power Lj. For given values of e e and Lj, while 
the value of (^Lv) pca k remains nearly constant, the fre- 
quency range, where the spectrum is flat, varies with the 
source size because of the changes in 7b r , 7max and the 
main emission mechanism (synchrotron or IC). It is em- 
phasized that the peak luminosity is chiefly governed by 
e e and Lj and is quite insensitive to the magnetic field 
strength and seed photons, which will only affect the fre- 
quency range of the flat spectrum. This means that if 
Lj is constrained by other inde pendent methods (e.g., 
lAllen et alj|2006t llto et al J 12051) . the observation of the 
peak luminosity will enable us to obtain information on 
the highly unknown acceleration efficiency e e . 

4. COMPARISON WITH EMISSIONS FROM THE LOBE 

Despite our focus on the non-thermal emission from 
the shell, emission from the lobe (or, cocoon) is an an- 



7 






1 A60 
1 U 




10 8B 




10°° 






10 57 


CM 




1 056 
1 U 




10 55 




1 D54 




1 U 




10 60 




10 58 








10 58 


(M 






10 57 





1 ' 1 ' 


I 1 


1 J 


r 








r 




t "s. 


-| 






I A 


-, 




L 45 =l 








R,= l 








i,i, 


1 





10 56 r 
10 55 i- 





Fig. 2. — Cooling and acceleration time scales (left panels) and energy distributions of non-thermal electrons (right panels) for sources 
with the jet power of Lj = 10 45 ergs s -1 . The top, middle and bottom panels are shown for the source sizes of R = 1, 10 and 100 kpc, 
respectively. The thick lines in the left panels give the total cooling time scale (thick solid line), the acceleration time scale (thick dashed 
line) and the source age (thick dot-dashed line) whereas the thin lines are contributions to the total cooling time scale from various processes: 
adiabatic losses (dot- dot- dashed line), synchrotron emissions (thin solid line), and IC scatterings of UV disc photons (long- short- dashed 
line), IR torus photons (thin dot-dashed line), NIR host galaxy photons (dotted line), lobe photons (thin short-dashed line) and CMB (thin 
long-dashed line) photons. 



other important ingredient which arises as a consequence 
of interaction of j et wi th ambient medium. As shown 
bv iStawarz et all (|2008l ). lobes can produce very bright 
emissions in broadband. Therefore, in order to consider 
the application of our model to the observation of radio 
sources, it is essential to investigate the relative signif- 
icance of the two emissions, since the lobe component 
may hamper the detection. In this section, we evaluate 
the emissions from the lobe and show the quantitative 



comparison between the spectra of the lobe and shell. 

4.1. Model for Emissions from the Lobe 

We basically follow the same procedure employed for 
the shell in calculating the energy distribution of the elec- 
trons and the resulting emissions. Since we are consider- 
ing the electrons residing in the lobe, the energy stored 
in the cocoon, E c = P c V c /{"f c — 1), is the energy budget 
for the emissions. From the dynamical model described 
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in ij2.ll the total internal energy stored in the cocoon is 
given by 



E c — fiobcLjt, 



(25) 



where / lobc = (5- a)(7 -2a)/[2a 2 + (1 -18%)a + 63%- 
28]. For 7e = 4/3, % = 5/3 and a = 1.5, /i obe = 7/13 ~ 
0.5 is obtained. As can be seen from the above equation 
and Eq. ([5]), E c is larger than E s by a factor of ~ 5. 

The energy distribution of the electrons in the lobe is 
determined based on Eq. (0 as in the case of the shell. 
As described in by evaluating the electron injection 
rate, Q(7 e ), and the cooling rate, 7cooi(7e), the energy 
distribution of the electrons is obtained by putting these 
quantities in Eq. (fl9| which corresponds to the solution 



of Eq. ([7]). In the present study, we assume that the 
electrons injected into the lobe have a power-law energy 
distribution given as 

Q(je) = K\ ohc J~ Pl ° b °f0I 7min,lobc < 7e < 7max,lobc- 

While the employed values of the power-law index and 
minimum Lorentz factor are the same as those adopted 
for the shell (piobo — 2 and 7min,lobe = 1), the max- 
imum Lorentz factor is fixed as 7max,iobe = 10 . We 
will comment on the assumed value of 7max,iobe later 
in $4.'2{ The normalisation factor, itjobc, is deter- 
mined from the assumption that a fraction, e ej i bc, of 
the energy deposited in the lobe is carried by non- 
thermal electrons: jT™> a *. lobe ^ _ l^m e c 2 Q( , y e )d'y e = 



9 



w 

00 
1) 



w 

OD 
S-. 

<D 

a 




9 12 15 18 21 24 27 
log[i/(Hz)] 




9 12 15 18 21 24 27 
log[i/(Hz)] 



Q43 
Q42 
Q41 
Q40 
Q39 
Q38 
Q37 



I 1 1 I 

R^lOO 



T 



I 1 1 I 1 1 I 1 1 I 1 1 



r 



i . . i . / 1 




'// \V\1 
• ./ 1 i i 1 1 



9 12 15 18 21 24 27 
log[u(Hz)] 



-J 




9 12 15 18 21 24 27 
log[u(Hz)] 



10 43 



- w 



w 
tafl 
s- 



-J 




9 12 15 18 21 24 27 
logf^Hz)] 



m 
<u 



10 43 
10 4S 
10 41 
10 40 
10 39 
10 38 
10 37 



'I 1 1 I 1 1 I 1 1 I 1 1 I 1 1 I 1 1 I 1 1 




9 12 15 18 21 24 27 
log[i/(Hz)] 



Fig. 4. — Broadband emission produced within the shell of sources with the jet powers of Lj = 10 45 ergs s 1 (left panels) and Lj = 
10 47 ergs s _1 (right panels). The top, middle and bottom panels are displayed for the source sizes of R = 1, 10, and 100 kpc, respectively. 
The various lines show the contributions from the synchrotron emissions (thin solid line) and IC scatterings of UV disc photons (long- 
short- dashed line), IR torus photons (dot-dashed line), NIR host galaxy photons (dotted line), CMB photons (long-dashed line) and lobe 
photons (short-dashed line). The thick solid line is the sum of these emissions. 



e e dE c /dt = /iobe£e,iobc£j- A rough estimation of 
the normalization factor is obtained as -fTiobe ~ 
0.5e ej i o bo£j/[m e c 2 ln(7 maX!lo bo)], where we used /i obG ~ 
0.5. Regarding the cooling rate of the electrons, we 
take into account the adiabatic losses as well as radia- 
tive losses due to the synchrotron and IC emissions. The 
adiabatic cooling rate is evaluated from Eq. (jTU]) and 
coincides with that of the shell. The synchrotron cool- 
ing rate is determined from Eq. under the assump- 
tion that a fraction ea, of the energy E c is carried by 
the magnetic fields. The corresponding energy density 



of magnetic field is given as 



e B E c 

V c . 



, , , , ,, i 1/3 j-2/3 n - 11/6 

4.0 x 10 M e B ,o. Wi L 45 R i er g cm 



where £b,o.i = cb/0.1. As for the IC cooling rate, we 
take into account all seed photon fields which are con- 
sidered in the shell except for the synchrotron photons 
from the lobe. Considering the IC of lobe photons, al- 
though a particular spectrum was assumed in the case 
of the shell (see Ej2.2[) . to be self-consistent, the distri- 
bution of the seed photons should be determined from 
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the calculated synchrotron emissions (synchrotron self- 
Compton). However, the method given in ^2 .21 used to 
solve Eq. ([7]) cannot be applied when a seed photon field 
has dependence on the electron distribution. Therefore, 
since the cooling rate due the synchrotron self-Compton 
is modest at best in any case considered in the present 
study, we neglected its effect merely for simplicity. As in 
the case of the shell, we assume that the photons from 
the disc, torus, host galaxy and CMB are monochromatic 
with frequencies given in £ 32.21 On the other hand, while 
same value is employed for CMB, the energy densities 
of the disc, torus and host galaxy photons are taken to 
be larger than those in the shell by a factor of 3, since 
the lobe is located closer to the emitting sources. The 
cooling rate due the IC of these photons are evaluated 
from Eq. (|12[) based on the above mentioned spectrum 
and energy densities. 

The emission spectrum of the lobe is calculated by fol- 
lowing the procedure described in £12.31 From the ob- 
tained energy distribution of the non-thermal electrons, 
spectra of the synchrotron and IC emission are evaluated 
from Eqs. (|20|) and (f2~T|) . respectively. The magnetic field 
and seed photon fields considered above is used for the 
calculation. In addition, although we did not consider 
its contribution on the cooling rate, we also evaluate the 
IC of the photons from the lobe based on the calculated 
synchrotron spectrum. 

4.2. Lobe vs Shell 

For an illustrative purpose, here we examine the emis- 
sions from the lobe by assuming e e! i bc = 1, an extreme 
case where all energy stored in the cocoon is converted 
to that of the non-thermal electrons, and £b = 0.1, a 
magnetic field strength factor of few below the equipar- 
tition val ue. These values are identical to those em- 
ployed in IStawarz et all ((2008). In Fig. [S] we show 
the obtained photon fluxes, vF v , from the lobe. In or- 
der to illustrate the comparison between the emissions 
from the lobe and the shell, photon fluxes from the shell 
are also shown in the figure. Same set of parameters 
which was previously assumed (Table. [T]) is employed 
for the emissions from the shell. It is assumed that 
the jet powers are Lj = 10 45 ergs s _1 {left panels) and 
Lj = 10 47 ergs s _1 (right panels) and the source is lo- 
cated at a distance of D = 100 Mpc. The top, mid- 
dle and bottom panels of the figure correspond to the 
source sizes of R = 1 kpc, 10 kpc and 100 kpc, re- 
spectively. We also plot the sensitivities of the Fermi 
satellite (http://www-glast.stanford.edu/) and HESS 
(http://www.mpi-hd.mpg.de/hfm/HESS/) and MAGIC 
(http://magic.mppmu.mpg.de/) telescopes. The dashed 
and solid lines display the photon fluxes from the lobe 
and shell, respectively. In addition to the total photon 
flux (thick black line), the contributions from the syn- 
chrotron emissions (thin black line) and the IC scatter- 
ings of UV disc photons (blue line), IR torus photons (red 
line), NIR host galaxy photons (light blue line), CMB 
photons (green line) and lobe synchrotron photons (pur- 
ple line) are also shown. 

As is seen in the figure, IC emissions from the lobe 
tend to be brighter when the source size is smaller as in 
the case of the shell. On the other hand, however, the 
emissions are synchrotron-dominated even for compact 
sources (R ~ 1 — 10 kpc) since the magnetic field strength 



is larger than that of the shell. As in the case of the 
shell, the luminosity scales linearly with e ei i o b and L } 
(see ^2.3|) . A rough estimate to the peak luminosity is 
obtained as 

(^WkJobc ~ 2.0 x 10 43 e e! i ob c^45 ergs s^ 1 . (26) 

It should be noted that the photon fluxes displayed in 
the figure correspond to an upper limit since e e j bo = 1 
is assumed. 

From Eqs. (|24|) and (|26|) . the ratio between the peak 
luminosities of the emissions from the lobe and shell is 
obtained as ~ 5e e jobeAe- Hence, the contrast between 
the two emissions is determined by the parameters e ej i bo 
and e e . If e e! i bc is much larger than e e , the overall spec- 
tra is dominated by the lobe. Although there is little 
constraint on the values of e ej i bo and e e , it is expected 
that the ratio e e ,iobc/ee is indeed large in most of the 
radio sources as is examined here due to the fact that, 
while large number of radio observations identified non- 
thermal emissions from the lobe, no clear evidence of 
the shell emissions have been reported so far at radio. 
It is emphasized, however, that, even in the case where 
£e,iobc is significantly larger than e e , there are frequency 
ranges in which the emissions from the shell can over- 
whelm those from the lobe as is seen in Fig. [5] This is 
due to the difference between the maximum energy of the 
electrons in the lobe and shell. Since 7 max ~ 10 s largely 
exceed 7 ma x,iobc = 10 5 , the synchrotron and IC emissions 
from the shell extend up to much higher frequencies than 
those from the lobe, and, as a result, the shell can dom- 
inate the emission spectra at the frequencies above the 
cut-off frequencies of the synchrotron and IC emissions 
from the lobe. 

Regarding the synchrotron emission, although the high 
frequency end of the emission is overwhelmed by the low 
frequency tail of the IC emissions from the lobe, emis- 
sions from the shell can become dominant at frequen- 
cies above the cut-off frequency of synchrotron emis- 
sion from the lobe which is roughly given as hv ~ 

6 x W-\ lm ^ lobe /10 5 )^l^Lli 3 R-^ eV when 

the source size is large (R > J?ic/ S yn)- F° r example, 
emissions from the shell can be observed at frequencies 
from IR to optical without strong contamination from 
the lobe for sources with size of R ~ 100 kpc (see the 
bottom panels of Fig. [5]). Hence, deep observations of 
radio galaxies at these frequencies may lead to discovery 
of the shell emissions. It is noted that a broader range 
of frequencies could be detectable if the ratio e e ,iobeAe 
is smaller. For compact sources (R < -Ric/syn)j however, 
synchrotron emissions are likely to be completely over- 
whelmed by the lobe components, since the synchrotron 
emissions from the shell is strongly suppressed (see the 
top panels of Fig. [5]). 

On the other hand, IC emissions from the shell are not 
subject to the contamination from the lobe irrespective 
to the size of the source at frequencies above the cut-off 
frequency of IC emissions from the lobe which is roughly 
given as hv ~ 7max,iob c m e c 2 ~ 50(7 maXj i ob c/10 5 ) GeV, 
since the lobe emissions cannot extend above the fre- 
quency. It is also noted that compact sources (R < 
.Rjc/syn) are favored for detection, since the luminosity 
is higher. Contrary to the above argument, although the 
same value was assumed for the maximum Lorentz fac- 
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tor (7 m ax,iobc = 10 5 L iStawarz et all (|2008T ) argued that 
the lobe can produce prominent emissions up to ~ TeV 
gamma-ray. It should be noted, however, that the emis- 
sions above the frequency hv ~ lOGeV in their study is 
considerably overestimated since the KN effect was not 
taken into account. Hence, due to the absence of lobe 
emissions, our study demonstrates that the energy range 
of hv > 10 GeV is the most promising window for the 
detection of the shell emissions. This can be confirmed 
in Fig. [5] indeed. For example, while the detection of 
~ GeV gamma-rays by the Fermi telescope may be ham- 
pered by the emissions from the lobe, ~ TcV gamma-rays 
are accessible to the MAGIC, HESS and, although not 
displayed in the figure, also VERITAS gamma-ray tele- 
scopes for the most powerful source with the jet power of 
Lj = 10 47 ergs s~ 4 located at D = 100 Mpc. It is worth 
noting that the maximum Lorentz factor 7 max of non- 
thermal electrons can be constrained by the detection of 
high energy cut-off at hv ~ 50(7 max /10 8 ) TeV in the 
IC emissions. Hence, the observation of the high energy 
gamma-rays can provide us with information not only on 
e e as mentioned in §3.21 but also on the gyro-factor £. 

Lastly, let us comment on the employed values of j max 
and 7max,iobc- In the above discussions, we have shown 
that, even in the case of e e ,iobc 3> e e , the high frequency 
part of synchrotron (~ IR — optical) and IC emissions 
(~ 10 GeV — 10 TeV) from the shell can overwhelm the 
lobe emissions. The above argument will be valid as long 
as (i) 7 max > 10 7 and (ii) 7 max ,iobo < 10 5 are satisfied. 
The condition (i) holds if the acceleration at the shell 
takes place in nearly Bohm limit (£ ~ 1) as shown in 
Figs. [5] and [3] As mentioned in §3J since the property of 
the shock considered here is similar to those of the SNRs 
in which electro n acceleration at nearly Bohm diffusion 
is inferred (e.g.. | Ya,mazaki et aLlfeOO^ iStage et al.l 120061 : 
iTanaka et al.l 120081 ) . it is natural to expect that this is 
also the case for the shell. In fact, £ ~ 1 inferred in 
the shell associated with radio galaxy Centaur us A from 
the X-ray observations (|Croston et al.l [2009). Hence, 
we expect that condition (i) is satisfied for the sources 
we focus here. Regarding the value of 7 roax ,lobe) multi- 
wavelength studies of powerful FRII radio galaxies show 
that maximum Lorentz factor of electrons residing in the 
lobe does not extend well beyond ~ 10 5 ((Stawarz et al.l 
120071: IGodfrev et al.ll2009t lYaii et alJl2010h Fl. Although 
the detail of the acceleration process is highly uncer- 
tain, if this feature is also common in compact sources 
(R ~ 1 — 10 kpc), condition (ii) is likely to be satisfied 
irrespective to the source size. It is wor th to note that 
the recent study bv lOstorero et al.l (|2010f ) has shown that 
the broadband spectra in some compact radio sources can 
indeed be fitted under the assumption of 7max.iobc ~ 10 5 . 

5. SUMMARY 

We have explored the temporal evolution of the emis- 
sions by accelerated electrons in the shocked shell pro- 



6 For "low-power" FRII radio sources with hotspots showing 
synchrotron emissions above optical freque ncies, however, elec- 
trons may be accelera ted to higher energies (Brunctti ct al. 2003; 
IHardcastle et al.|[200l 1. 



duced by AGN jets. Focusing on the powerful sources 
which host luminous quasars, we have calculated the 
spectra of the synchrotron emission as well as the IC 
scatterings of various photons that will be relevant in 
this context. We have used a simple analytic model that 
describes the dynamics of the expanding shell and esti- 
mated the energy distribution of non-thermal electrons 
based on this model, taking properly into account both 
the adiabatic and radiative coolings. Below we summa- 
rize our main findings in this study. 

1. When the source is small (R < Ric/ syn ~ 
27L^ 46 BZ^pc), the dominant radiative process is the 
IC scattering of IR photons emitted from the dust torus. 
For larger sources, on the other hand, the synchrotron 
emissions dominate over the IC emissions, since the en- 
ergy density of photons becomes smaller than that of 
magnetic fields (Ub > U p h cx R~ 2 ). Through the entire 
evolution, the spectrum is rather broad and flat, and the 
peak luminosity is approximately given by {vL v ) vca ^ ~ 
4.0 x 10 40 e_2^45 ergs s _1 , since it is roughly equal to the 
energy injection rate, which is in turn determined by the 
jet power Lj and acceleration efficiency e e . 

2. The broadband spectra extend from radio up to 
~ 10 TeV gamma-ray energies for a wide range of source 
size (R ~ 1 — 100 kpc) and jet power (Lj ~ 10 45 — 
10 47 ergs s" 1 ). By comparing the emissions with those 
from the lobe, we find that the synchrotron emissions at 
IR and optical frequencies can be observed without be- 
ing hampered by the lobe emissions for extended sources 
(R > i?ic/syn)i while the IC emissions at hv > 10 GeV 
can be observed with the absence of contamination from 
the lobe irrespective of the source size. In particular, 
it is predicted that, for most powerful nearby sources 
(Lj - 10 47 ergss" 1 , D < 100 Mpc), - TeV gamma- 
rays produced via the IC emissions can be detected by 
the modern Cherenkov telescopes such as MAGIC, HESS 
and VERITAS. 
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